Self-consistent spin waves in magnetized BEC 
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We obtain equations of quantum hydrodynamic (QHD) for magnetized spin-1 neutral Bose-Einstein 
condensate (BEC). System of QHD equations contains the equation of magnetic moment evolution (an 
analog of the Bloch equation). We account spin-spin interaction along with the short range interaction. 
We consider self-consistent field approximation of QHD equations. Starting from QHD equation we derive 
the Gross-Pitaevskii equation for magnetized BEC. We show that Gross-Pitaevskii equation exists under 
condition that the magnetic moment direction is not change. Using obtained QHD equations we study the 
dispersion of collective excitation. As in electrically polarized BEC [P. A. Andreev, L. S. Kuz'menkov, 
arXiv: 1201.2440], in the magnetically polarized BEC there is second wave mode (polarization mode or 
spin wave), in addition to the Bogoliubov's mode. Second wave solution appears due to the magnetic 
moment evolution. The influence of magnetization on dispersion of Bogoliubov's mode is found. We 
found strong difference of dispersion properties of waves in magnetized BEC from electrically polarized 
BEC. 
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I. INTRODUCTION 

Special attention at studying Bose-Einstein condensation 
(BEC) in ultracold atomic vapors gives to polarized BEC 
El- uM- The term a polarized BEC is used in literature as 
for electrically HI, HI as for magnetically H3l, [QJ] 
polarized BEC. In the recent years, an attention increases 
to the BEC of atoms with magnetic polarization because of 
the BEC of 52 Cr was experimentally realized 11711 . These 
atoms have a large magnetic moment 6/i#, where \i B is 
the Bohr magneton. There are a lot of attempts to realize 
electrically polarized BEC 03], lfl6Tl . for this aim two 
nuclear molecules are usually used. Attention to electri- 
cally polarized BEC is particularly caused by the fact that 
molecules have electrical dipole moment of order of 1 De- 
bye and this means that the force of dipole-dipole inter- 
action of molecules on four orders more than the force of 
spin-spin interaction of such atom as 52 Cr. 

Nevertheless, in the chromium atoms the quantity of 
spin-spin interaction is comparable with the value of the 
short range interaction. Thereby, a spin-spin interaction 
makes anisotropy. Consequently, spin- spin interaction has 
essential contribution in chromium atoms static and dy- 
namic properties. 

Usually, for polarized BEC studying a generalized 
Gross-Pitaevskii (GP) equation is used [l]] - ifllll . A non- 
linear Schrodinger equation is usually used for quantum 
gases studying. One is also used for dipoles BEC descrip- 
tion. In the last case, in GP equation a new term is added. 
However, this approximation includes interaction among 
parallel dipoles, and does not account spatial and tempo- 
ral dipoles evolution. Different approximations based on 
scattering process were analyzed in Ref. ifl^l . Particu- 
larly, where were considered condition there can be used 
first Born approximation and presented generalization of 
GP equation for the scattering of polarized atoms beyond 
first Born approximation. In Ref. [12] authors devel- 
oped method accounted spatially inhomogeneous distribu- 



tion of dipoles directions. This method includes the poten- 
tial of interaction unparallel dipoles in GP equation, and the 
second equation determining dipoles moment of a volume 
unit. In paper lfl3ll were suggested a three fluid hydrody- 
namic model for description of a spin-1 Bose gas at finite 
temperatures, three sounds are evaluated there in a mean- 
field approximation. The method of many-particle quan- 
tum hydrodynamic for electrically polarized BEC was de- 
veloped in Ref. |fl4ll . This approximation allows describe 
as spatial as temporal evolution of electric dipole moment. 
For magnetized BEC studying in this paper we use the 
method of quantum hydrodynamics (QHD). This method 
was developed in the recent years for various physical sys- 
tems. Namely, it was quantum plasma 11201 . Il2ll l.plasma 
of particles with the own magnetic moment l22ll- I128ll . rel- 
ativistic quantum plasma 12911 . ultracold Bose and Fermi 
gases with nonlocal interaction II30|] and three-particle in- 
teraction at nonzero temperatures 13111 . quantum particles 
with electrical polarization ll32ll . particularly being in the 
BEC state lfl4ll . D3I1 . the graphene electrons 13411 and 
BEC of graphene excitons [35]. Influence of magnetic 
moment evolution on the collective excitation dispersion 
for quantum plasma were considered in Ref. [361], where 
were found that temporal evolution of magnetic moment 
leads to existence of new branches in dispersion depen- 
dence. More detailed analysis of spin waves in quantum 
plasma and methods of these waves generation were con- 
sidered in Ref. JH. 

In this paper we develop a self-consistent theory for a 
magnetic moment evolution of the BEC of neutral spin-1 
atoms. We derive the QHD equations: continuity equation, 
momentum balance equation (Euler equation) and mag- 
netic moment evolution equation. The last one describe 
the temporal evolution of the BEC magnetization and it's 
influence on BEC properties. In Ref. ifbUl were shown 
that temporal evolution of collective electric polarization 
in BEC of particles having electric dipole moment leads to 
existence of new type waves in BEC - polarization waves, 
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where supposed that BEC located in external uniform elec- 
tric field. Dispersion dependence of collective excitations 
shows no dependence on direction of wave propagation and 
exhibits comparably stable excitation spectrum. Here, we 
study the linear dispersion in three dimensional magnetized 
BEC located in external magnetic field. We find the contri- 
bution of equilibrium magnetization in dispersion depen- 
dence of Bogoliubov's mode and obtain that magnetization 
leads to destabilization of Bogoliubov's mode. The last fact 
coincide to the other authors results f5D- J8fl. We show ex- 
istence of magnetic moment waves (spin waves). The spin 
waves propagate along direction of external magnetic field. 

Our paper is organized as follows. In Sec. II we de- 
scribe the used model and present the QHD equations for 
magnetized BEC. One of the equations is the magnetic 
moment evolution equation, this equation describe evolu- 
tion of magnetic moment in space and time due to the in- 
teraction particles magnetic moments each other and with 
the external field. Corresponding GP equation is also pre- 
sented in Sec. II. Condition for GP equation derivation is 
described. In Sec. Ill we presented the method of QHD 
equations solution and the method of obtaining of the dis- 
persion dependence. In Sec. IV the dispersion spectrum 
properties of the collective excitations are described. In 
Sec. V brief summary of obtained results is presented. 



II. MODEL 

We do not consider here the details of the QHD equation 
derivation, but we present and use in the paper the QHD 
equation for the neutral Bose particles with the own mag- 
netic moment being in the BEC state. All details for equa- 
tion derivation can be found in Ref.s lEoll . 1221. OOh . 
1 32j]. Here we present the Schrodinger equation used for 
QHD equation derivation 



ihd t ip s (R,t) 




Vl _ rsl. £ a U a 

2 B i[ext] 



j^Gtf^s^U) (R,t), (1) 



where we include short range and spin-spin interactions, 
and action of external magnetic field on particles spin. In 
the Schrodinger equation fl]) we use following designa- 
tions: 7i is the gyromagnetic ratio, p 4 = —ihdi is the op- 
erator of momentum, B^e^-is the external magnetic field 
acting on i-th particle, Uij present the short range interac- 
tion (SRI), the Green function of the spin-spin interaction 
(SSI) has form 



G?f = d° ! d?(l/r ij ) + 4Tr5 a(3 6(r ij ). 
For spin matrixes the commutation relations are 



(2) 



Thereby, we consider spin-1 Bose particles we present 
here the evident form of the spin matrixes sf for particles 
with spin equal to 1 : 
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At derivation of system of QHD equations we need to 
write an evident form of Hamiltonian of dipole-dipole in- 
teraction. Usually using expressions for a Hamiltonian of 
dipole-dipole interaction are the same as for electric as for 
magnetic dipoles: 

Hdd = z , (3) 



but in this paper for SSI we used ©. In Ref. fl 1 411 were 
shown that for electric dipoles interaction we need to use 
the first term in formula (O. 
There is well-known identity 



d a d p - 
r 



_ q a ji 12 A 

i ^_L_LL + Zls a0 d(r), (4) 



so we can see the difference between usually using Hamil- 
tonian and one's used in this paper (0. The corrections 
of our choice followed from the fact that the equations ob- 
tained in the paper coincide to the Maxwell equations. 

It has been shown by Breit [37]- li40h that a Hamilto- 
nian for spin-spin interaction, and, as a consequence, for 
the interaction of magnetic moments contains a term that is 
proportional to Dirac <5-function 5(ri—r 2 )d"d,2 along with 
©. The coefficient of the (^-function has been refined later 
I22II by using idea that the Hamiltonian is in accord with 
Maxwell's free equations, such as divB = 0. The resultant 
expression for the spin-spin interaction Hamiltonian is: 



(5) 



Formula (01 is obtained at substitution of relation d4]l in to 
the formula 

The system of QHD equations for magnetized BEC con- 
sists of: 
continuity equation 

d t n{r, t) + d a {n{r, t)v a (r, t)) = 0, (6) 

where n(r, t) is the concentration of particles, v(r, t) is 
the velocity field. In this equation describing the whole 
particle number invariance, and evolution of concentration 
owing to particles displacement; 

momentum balance equation (Euler equation) 

mn(r, t) (d t + vV)/ (r, t) + d p p a0 (r, t) 
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4m 4m \ n(r, t) 



Tn(r,t)d a n(r,t) + -T 2 d a An 2 (r, t) 



+ M^(r,t)d a B^(r, t) 



(7) 



where V is the the gradient operator, A is the Laplace op- 
erator, M(r, t) is the density of magnetic moment (mag- 
netization), B(r, t) is the magnetic field, p a "(r,t) is the 
kinetic pressure describing the temperature effects and can 
be vanished for the BEC, 



and 
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T 
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dr{r) 



3 dU(r) 



dr 



dr{r) 



5 8U(r) 



dr 



(8) 



(9) 



describing the SRI influence on particles dynamic. The 
momentum balance equation for the BEC with SRI up to 
TOIR approximation were obtained in Ref. 13011 in the ab- 
sence of the magnetization contribution. The second term 
in the right-hand side of equation (Q is the term appear- 
ing in the TOIR approximation. This term is an example 
of non-local SRI, analogous terms were obtained in Ref.s 
111, JH; 

balance equation of magnetic moment (Bloch equation) 

d t M a (r,t) + V p J?f (r,t) = -e^M^tr, t)B~<(r, t), 

n 

(10) 

where a tensor of spin current J^f arises. Vanishing by 
thermal motion we have J^f = M a v f3 . This equation de- 
scribes the dynamic of magnetization and influence of ex- 
ternal magnetic field and inter-particle interaction on mag- 
netization evolution. For the first time the many-particle 
QHD equations for spinning particles were derived in Ref. 
1 220 for the case spinning charged Fermi particles. It was 
made in 2001. 

During the derivation of QHD equation we do not con- 
sider the scattering problem and do not use conception of 
scattering for interpretation of the SRI and SSI. However, 
for simplicity of comparison of our results with the results 
of other authors we reduce connection of T and T 2 with 
the scattering amplitude [13011 . 

The first order of the interaction radius interaction 
constant for dilute gases has the form T = — g = 
—4irh 2 a/m, where g is the interaction constant usually 
used in GP equation, a is the scattering length ll30l l43ll . 
The value Y 2 may be expressed approximately as T 2 = 
—8a 2 T/8 [30], where 9 is a constant positive value about 
1, which depends on the interatomic interaction potential. 
Finally, T 2 takes the form T 2 = — Tr9h 2 a 3 /2m. 

We consider the three-dimensional (3D) system of Bose 
particles. It allows us to introduce in equations (0 and 



([Tol l the magnetic field caused by particles magnetic mo- 
ments. Spin-spin interaction in these equations at deriva- 
tion appears in the integral form (see for example [22]). 
In self-consistent approximation in 3D particles system we 
can represent interaction via magnetic field which has form 

B°(r,t) = J dr'G a0 {\ r-r' \)M (r',t), 

which satisfies to the Maxwell equation: 

curlB(r,t) = 4ncurlM(r,t). (11) 

From equations (l6ll-(fT0ll we can see that short range in- 
teraction, the terms proportional to T and T 2 , gives contri- 
bution in the momentum balance equation only. 

In the absence of magnetic moment direction evolution 
M a (r,t) = n(r, t)/j, a (r, t), // Q (r, t) = const (suppos- 
ing that magnetization change because of the concentration 
changing only), we can derive the GP equation for spinning 
neutral particles in the BEC state. We notice that for the 
fi a (r,t) = const the magnetic moment evolution equation 
does not exist (It does not need). In this case the magnetic 
moment evolution reduces to the concentration evolution 
M(r, t) ~ n(r, t). Particles dynamic is described by two 
equations: the continuity equation © and simplification of 
the Euler equation ©. 

The method of GP equation derivation from QHD equa- 
tion is described in Ref.s Q and JH. Resulting GP 
equation, in the absence of terms appeared up to TOIR ap- 
proximation, has form 



ihd t ®(r,t) 



2m 



V 2 + g | $(r,t) 



+ l4ti J dr'G^(r,r')|$(r',t)| 2 ^(r,t), (12) 

where G afi the Green function of SSI ©, using evident 
form of the Green function we can rewrite GP equation 

iHdMr,t) = l-^-V 2 + g\ $(r,t) | 2 



■"> / -, 1 1 3 cos \ i q 

+H 2 / dr'— ^|$( r ',t) 12 



r — r 



'13 



ftimr,t)\ 2 U(r,t) 



(13) 



This equation has difference from usually used for mag- 
netized BEC H, Jl, U, OD. This difference appear 
because of second term in ® proportional to delta func- 
tion. The last term in equation (fT3l) appear because of the 
second term in Green function of SSI 

At low temperature in the external magnetic field the 
atoms magnetic moments directed along the external field, 
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because of absence disordered factors as a thermal motion. 
However, the perturbation of dipoles direction might prop- 
agate there. These perturbations can exhibits as waves. In 
this case we have deal with the spin waves. It is very spread 
phenomenon in the system of particles with magnetic mo- 
ment. For example we can admit that spin waves there are 
in ferromagnetic J4J], but there magnetic moments bound 
with the unit of crystalline lattice. In gases atoms can move 
through the system. Thus, in gases there is one more mech- 
anism of magnetization changing, compare with ferromag- 
netic. On the other hand, the Bose-Einstein condensation 
of spin wave quanta has been also studied 14511 . 

Different generalization of GP equation has been used in 
literature. One of them is the spinor nonlinear Schrodinger 
equation 14611 . This equation is equivalent to the three hy- 
drodynamic equations ll22ll . these are continuity equation, 
Euler equation and magnetic moment balance equation. 
Whereas, usual scalar GP equation equivalents to continu- 
ity and Euler equations only. Thus, the spinor generaliza- 
tion of GP equation also described the magnetic moment 
evolution. Consequently, this equation contains informa- 
tion about spin waves considered in our paper. At deriva- 
tion of the hydrodynamic equations from the spinor nonlin- 
ear Schrodinger equation the SSI appears in self-consistent 
field approximation, without two-particle quantum correla- 
tion. In fact, quantum spin-spin correlations might be in- 
cluded in GP equation if they have a simple enough form. 

In the case of spin-0 particles or if we do not account 
magnetic moments motion there is the collective excitation 
in the BEC, it is the Bogoliubov mode. In Ref. [14] were 
shown that electric dipoles evolution leads to the existence 
of new wave mode, it is the wave of electrical polariza- 
tion. In this paper, we show the existence of the spin waves 
(magnetic moment waves) in magnetized BEC. Below we 
consider wave dispersion calculations based on above de- 
scribed model. 

IH. METHOD OF SOLUTION 

We can analyze the linear dynamics of elemental excita- 
tions in the magnetized 3D BEC using the QHD equations 
©, (0, ( TTOb and ( fTTT) . Let's assume the system is placed 
in an external magnetic field B = B e z . The values of 
equilibrium concentration n and equilibrium magnetiza- 
tion M ~ B for the system in an equilibrium state are 
constant and uniform and its velocity field v Q (r, t) values 
equal to zero. 

We consider the small perturbation of equilibrium state 
n = n Q + Sn, v a = + v a and M a = M " + 5M a . Sub- 
stituting these relations into system of equations ©, (0, 
( TTOb and ( fTTT) and neglecting nonlinear terms, we obtain a 
system of linear homogeneous equations in partial deriva- 
tives with constant coefficients. Passing to the following 
representation for small perturbations 5f 

5f = f(oj, k)exp(—iujt + ikr) 



one yields the homogeneous system of algebraic equations. 
The magnetic field strength is assumed to have a nonzero 
value. Expressing all the quantities entering the system of 
equations in terms of the magnetic field, we come to the 
equation A Q/3 (w,k) • B"(w,k) =0, where A Q/3 (w,k) is 
the dispersion matrix. The evident form of the dispersion 
matrix A Q ^(w,k) is presented in the Appendix. In this 
case, the dispersion equation is detA(uj, k) = 0. Solv- 
ing this equation with respect to uj 2 we obtain following 
results. 



IV. ELEMENTARY EXCITATIONS IN THE POLARIZED 

BEC 

The dispersion characteristic for EE in BEC can be ex- 
pressed in the form of 

c 2 =(^ + ^y-^(T + 4^ 2 , (14) 
\4m^ my m 

the equilibrium magnetization /x caused by external mag- 
netic field B and has form /i = x^o/ ra o> where % 
is the ratio between equilibrium magnetic susceptibility 
and magnetic permeability. In the case when all magnetic 
dipoles directed parallel to external field (in equilibrium 
state) we have fi = 7. From formula ( fT4l) we can see that 
for the repulsive SRI (T < 0) at long wave length limit 
magnetization lead to instability of the Bogoliubov's mode 
because the frequency square become negative: 

^ = T I -Ait^k 2 < 
m 

at condition 47t//q >| T |. 

We obtained three solutions for spin waves. They fre- 
quencies are 

uj = f2, 

w = 0(l±4vrx), (15) 

where Q = -fB /his the cyclotron frequency. It is the fre- 
quency of the single magnetic moment precession in the 
external uniform magnetic field. For the case when all 
magnetic dipoles directed parallel to external field fully 
magnetized BEC, which can be considered as ferromag- 
netic state, we have 

/0(1 + 8tt) \ 
uj = n(l±4n X ) = I ^ I . (16) 

Thus, we have only two solution in this case because one 
of solutions uj = Sl(l ± 47rx) equal to 0. 



V. CONCLUSION 



We researched dispersion properties of magnetized 3D 
BEC. We obtained the contribution of equilibrium magne- 
tization in the Bogoliubov's mode dispersion dependence. 
We found that time evolution of magnetic moments lead 
to existence of new type of waves in magnetized BEC, 
it is spin waves. We found that SRI has no influence of 
spin wave dispersion. At comparison dispersion properties 
magnetized BEC with the electrically polarized BEC we 
got large difference between them. 

We studied dispersion of collective excitations by means 
of QHD equations derived in this paper for magnetized 
BEC. 

We have found two types of branches of dispersion 
curve. These are Bogoliubov's mode and spin waves. 

In correspondence with the other authors we obtained 
that magnetization leads to destabilizing of Bogoliubov's 
mode. For attractive SRI at small wave vector we have 
to 2 = — n | T | jm ■ k 2 = —AttH 2 \ a j n /m 2 ■ k 2 , this 
spectrum is unstable. Addition of magnetization leads to 
instability growing to 2 = — (n /m)(| T | +47r / u 2 l )£; 2 . 
For repulsive SRI T = — | T | and in the absence 
of the magnetization we have stable sound wave to 2 = 
n | T | /m ■ k 2 , but magnetization disturbs stability 
to 2 = (n /m)(| T | —4irfil)k 2 . This spectrum remains 
stable in two cases. The first one while magnetization is 
small compare with the constant of SRI T. The second 
case corresponds to the large constant of SRI T which can 
be reached by the Feshbach resonance |47|]- 04 811 . 



In general case in magnetized BEC exist three spin 
waves with frequencies: uj = O and to = 0(1 ± 47rx). 
But, for fully magnetized ferromagnetic BEC 0(1 — 47rx) 
equal to 0, so we have only two solutions to = Vt and 
oj = 0(1 + 8tt). As any physical system of particles with 
magnetic moments, especially ferromagnets, the magne- 
tized BEC shows existence of the spin waves. In this paper 
we consider the macroscopic spin dynamic described by 
QHD, but there are a lot of well-known method's based on 
second quantization Il49ll - 15311 - These methods are close 
to the Bese-Habbard model and can be used together for 
studying of microscopic influence of magnetic moment dy- 
namic on quantum gases properties. 
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APPENDIX: EVIDENT FORM OF DISPERSION 
MATRIX 

We describe A a/3 (uj,k) by presenting of the evident 
form of each component of the matrix 

7Mo 1 



A xx (oj,k) = 4:7rtok z n - 



h oj 2 -^ 



A xy (to,k) = ikg + 4irik z n 



7.60 7^o 1 

h h UJ 2_ ilS 



A xz (u),k) = — iky + 'imk y noH E; 
A^(a;,k) = (A*>,k))*; 
A yy (w,k) = A xx (to,k); 
A yz (to,k) = %k x — Amk x nofMoS; 
A zx {co,k) = iky 



+47rn 



7Mo 1 



tok x + ik. 



7^o 



n C0 2_ts\ x 1 v n ) ' 



A zy (to,k) = -ik. 



-47rn -^T-^ tttt ( —iok„ — ik ' 



h oji-i^lK v x h 7' 



where 



A zz (co,k) = 0, 

-jJL k 2 



muj 2 - £fe 4 + (T + 47r^)n F - n Q T 2 k^ ' 
Quantity H contains whole information about SRI. 
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